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Measurement of the D*(2010)~^ meson width and the Z?*(2010)"'" — Z?° mass difference 
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Measurements of the mass difference Amo between the _D*(2010)^ and the D" mesons and of 
the natural hue width V of the transition Z)*(2010)''' -^ D^tt'^ are presented. The data were 
recorded with the BABAR detector at center-of-mass energies at and near the T{4S) resonance, and 
correspond to an integrated luminosity of approximately 477 fb~^. The D" meson is reconstructed 
in the decay modes Z?" — ^ K~-k^ and Z?" — ^ K~-k^tt~tv'^ . For the decay mode D" -^ K~tt'^ 
we obtain F = (83.4 ± 1.7 ± 1.5) keV and Amo = (145 425.6 ± 0.6 ± 1.8) keV, where the quoted 
uncertainties are statistical and systematic, respectively. For the D" — > K~tv'^'k~tt'^ mode we obtain 
F = (83.2 ± 1.5 ± 2.6) keV and Amo = (145 426.6 ± 0.5 ± 2.0) keV. The combined measurements 
yield F = (83.3 ± 1.3 ± 1.4) keV and Amo = (145 425.8 ± 0.5 ± 1.8) keV; the result for the width 
is approximately 12 times more precise than the current world average, while that for the mass 
difference is more precise by a factor of approximately 6. 

PACS numbers: 13.20.Fc, 13.25Ft, 14.40.Lb, 12.38.Gc, 12.38.Qk, 12.39.Ki, 12.39.Pn 



The line width of the D*{2010)+ (D*+) provides a win- 
dow into a nonperturbative regime of strong interaction 
physics where the charm quark is the heavier meson con- 
stituent [l|-l3|. The line width provides an experimental 
test of D meson spectroscopic models, and is related to 
the strong coupling of the D*Dtt system, go'DTv In the 
heavy-quark limit, which is not necessarily a good ap- 
proximation for the charm quark [J], this coupling can 
be related to the universal coupling of heavy mesons to 
a pion, g. Since the decay B* — >■ B-k is kincmatically for- 
bidden, it is not possible to measure the coupling gB*BTT 
directly. However, the D and B systems can be related 
through g, allowing the calculation of gB*B-n, which is 
needed for a model- independent extraction of \Vub\ la,@| 
and which forms one of the larger theoretical uncertain- 
ties for the determination of \Vub\ H- 

We study the D*"*" — >■ Z)°7r+ transition, using the 
D" -^ K~TT'^ and D" -^ K^tt^-k'tt^ decay modes, to 
extract values of the Z?*+ width F and the difference be- 



tween the D*^ and D° masses Amp. Values are reported 
in natural units and the use of charge conjugate reactions 
is implied throughout this paper. The most precise prior 
measurement of the width [8[, by the CLEO Collabo- 
ration, is F = (96 ± 4 ± 22) keV, where the uncertainties 
are statistical and systematic, respectively. In the present 
analysis, we use a data sample that is approximately 50 
times larger. This allows us to apply restrictive selection 
criteria to reduce background and to investigate sources 
of systematic uncertainty with high precision. 

To extract F, we fit the distribution of the mass differ- 
ence between the reconstructed -D*+ and the J5° masses, 
A?7i. The signal component is described with a P-wave 
relativistic Breit-Wigner (RBW) function convolved with 
a resolution function based on a Geant4 Monte Carlo 
(MC) simulation of the detector response [9[. The de- 
tector resolution (full-width at half-maximum) for A?7i 
is approximately 300 keV, which is larger than the mea- 
sured value of F [8| . The key to this analysis is that, while 



the full-width at half-maximum of the almost-Gaussian 
resolution is larger than the extracted width, the RBW 
tails allow us to use information far from the central sig- 
nal region to measure the width. Consequently, we pay 
particular attention to how the selection criteria affect 
the tails. 

This analysis is based on a data set corresponding 
to an integrated luminosity of approximately 477 fb~^ 
recorded at, and 40 MeV below, the T{4S) resonance. 
The data were collected with the BaBAR detector at the 
PEP-II asymmetric energy e+e" collider, located at the 
SLAC National Accelerator Laboratory The BaBAR de- 
tector is described in detail elsewhere [10| ; we summarize 
the relevant features below. The momenta of charged 
particles are measured with a combination of a cylin- 
drical drift chamber (DCH) and a 5-layer silicon vertex 
tracker (SVT), both operating within the 1.5 T mag- 
netic field of a superconducting solenoid. Information 
from a ring-imaging Cherenkov detector is combined with 
specific ionization (dE/dx) measurements from the SVT 
and DCH to identify charged kaon and pion candidates. 
Electrons arc identified, and photons measured, with a 
CsI(Tl) electromagnetic calorimeter. The return yoke of 
the superconducting coil is instrumented with tracking 
chambers for the identification of muons. 

We remove a large amount of combinatorial and B me- 
son decay background by requiring D*+ mesons produced 
in e"'"e^ — > cc reactions to exhibit an e+e~ center-of- 
mass-frame momentum greater than 3.6 GeV. The entire 
decay chain is fit using a kinematic fitter with geometric 
constraints at the production and decay vertex of the D° 
and the additional constraint that the D*'^ laboratory 
momentum point back to the interaction region of the 
event. The pion from £)*+ decay is referred to as the 
"slow pion" (denoted 7r+) because of the limited phase 
space available in the D*~^ decay. The selection crite- 
ria are chosen to provide a large signal-to-background 
ratio (S/B), in order to increase the sensitivity to the 
long signal (RBW) tails in the Am distribution; they are 
not optimized for statistical significance. The criteria 
are briefly mentioned here and presented in detail in the 
archival reference for this analysis |11| . The experimen- 
tal uncertainty in Am is dominated by the uncertainty 
on the measured momentum of the slow pion. We imple- 
ment criteria to select well-measured pions. We define 
our acceptance angle to exclude the very-forward region 
of the detector, where track momenta are not accurately 
reconstructed, as determined using an independent sam- 
ple of reconstructed Kg ^ 7r~7r+ decays. The Kg re- 
constructed mass is observed to vary as a function of the 
polar angle 9 of the K'^ momentum measured in the lab- 
oratory frame with respect to the electron beam axis. To 
remove contributions from the very-forward region of the 
detector we reject events with any Z?*+ daughter track 
for which coa9 > 0.89; this reduces the final samples by 
approximately 10%. 
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FIG. 1. (color online) Fits to data for the D" — ^ K~tv'^ and 
D'^ — >■ K^Tv'^Tv^Tv'^ decay modes. The total probability den- 
sity function (PDF) is shown as the solid curve, the convolved 
RBW-Gaussian signal as the dashed curve, and the back- 
ground as the dotted curve. The total PDF and signal com- 
ponent are indistinguishable in the peak region. Normalized 
residuals are defined as (A'"obscrvcd - A''prcdictcd) /\/A^prcdictcd. 



Our fitting procedure involves two steps. In the first 
step we use the sum of three Gaussians and a threshold- 
hke function [11| to model the finite detector resolu- 
tion associated with track reconstruction, by fitting the 
Am distribution for correctly reconstructed MC events. 
These simulated D*^ decays are generated with F = 
0.1 keV, so that the observed spread of the MC distri- 
bution can be attributed to event reconstruction. The 
threshold function describes 7r+ decays in flight to a /i, 
for which coordinates from both the tt and fi segments 
are used in track reconstruction. 

The second step uses the resolution shape parameters 
from the first step and convolves the Gaussian compo- 
nents with a RBW function to fit the measured Am dis- 
tribution in data. The RBW function is defined by 
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where Td'Dtt is the partial width to D'^irf, m is the 
D^TT'^ invariant mass, toq is the invariant mass at the 
pole, and rTotai(w) is the total D*"*" decay width. The 
partial width is defined by 
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Here i = \, T\^^ (p) = y/TTr^ is a Blatt-Wcisskopf 
form factor for a vector particle with radius parameter r 
and daughter momentum p, and the subscript zero de- 
notes a quantity measured at the mass pole rn n Il2l [13 1 . 



We use the value r = 1.6GcV~ from Ref. [1^. For 
the purpose of fitting the Am distribution, we obtain 
dr{Am.)/dAm from Eqs. ([T]) and ^ through the substi- 
tution m = m{D'^) + Am, where m{D^) is the nominal 
D" mass ^. 

As in the CLEG analysis [8[ , we approximate the total 
D*^ decay width FTotai(?Ti) w T 0*0-^(^1), ignoring the 
electromagnetic contribution from _D*+ — )■ Z?+7. This 
approximation has a negligible effect on the extracted 
values, as it appears only in the denominator of the RBW 
function. 

To allow for differences between MG simulation and 
data, the root-mean-square deviation of each Gaussian 
component of the resolution function is allowed to scale 
in the fit process by the common factor (1 + e). Events 
that contribute to the non-Gaussian component have a 
well- understood origin (tt^ decay in flight), which is accu- 
rately reproduced by MG simulation. In the fit to data, 
the threshold function has a fixed shape and relative frac- 
tion, and is not convolved with the RBW. The relative 
contribution of the threshold function is small (< 0.5% of 
the signal), and the results from fits to validation signal- 
MG samples are unbiased without convolving this term. 
We fit the Am distribution from the kinematic threshold 
to Am = 0.1665 GcV using a binned maximum likelihood 
fit and an interval width of 50 keV. 



In the initial fits to data, we observed a strong depen- 
dence of Atoq on the slow pion momentum. This depen- 
dence, which originates in the modeling of the magnetic 
field map and the material in the beam pipe and SVT, is 
not replicated in the simulation. Previous _BaZ3AR stud- 
ies have observed this effect, for example the measure- 
ment of the A^ mass [16|. In that analysis, the material 
model of the SVT was altered in an attempt to correct 
for the energy loss and the under-represented small-angle 
multiple scattering (due to nuclear Goulomb scattering). 
However, the momentum dependence could only be re- 
moved by adding a nonphysical amount of material to 
the SVT. The corrections used here avoid altering the 
underlying detector model and focus on adjusting track 
momenta following reconstruction. We use a sample of 
K'g — ?> 7r+7r~ events from Z)*+ — >■ D^tt'^ decays, where 
D" — > /ir^7r+7r~, and require that the Kg daughter pi- 
ons satisfy the same tracking criteria as the 7r+ candi- 
dates for the D° —>■ K~Tr+ and £)° — >■ i4r~7r+7r^7r+ sig- 
nal modes. The Kg decay vertex is required to lie inside 
the beam pipe and to be well separated from the I?° ver- 
tex. These selection criteria yield an extremely clean Kg 
sample (over 99.5% pure), which is used to determine 
fractional corrections to the overall magnetic field and to 
the energy losses in the beam pipe and, separately, in the 
SVT. We determine the best set of correction parameters 
by minimizing the difference between the tt+tt^ invari- 
ant mass and the current world average for the K^ mass 
(497.614 ± 0.024 MeV) ^ in 20 intervals of laboratory 
momentum in the range 0.0 to 2.0 GeV. The parameters 
increase the magnitude of the magnetic field by 4.5 Gauss 
and increase the energy loss in the beam pipe and SVT 
i% and 5.9%, respectively 



111. The momentum- 
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dependence of Am,Q in the preliminary results was mostly 
due to the slow pion. However, the correction is applied 
to all £)*+ daughter tracks. All fits to data described in 
this analysis are performed using masses and Am values 
calculated using corrected momenta. Simulated events 
do not require correction because the same field and ma- 
terial models used to propagate tracks are used for their 
reconstruction. 

Figure[l]presents the results of the fits to data for both 
the 1?° -> K~Tr^ and -D" — >■ K~TT~^TT~Tr~^ decay modes. 
The normalized residuals show good agreement between 
the data and our fits. Table U summarizes the results of 
the fits to data for both 1?° decay modes. The table also 
shows S/B at the peak and in the high Am tail of each 
distribution. 

We estimate systematic uncertainties related to a va- 
riety of sources. The data are divided into disjoint sub- 
sets corresponding to intervals of D*+ laboratory mo- 
mentum, D*~^ laboratory azimuthal angle 0, and recon- 
structed D'^ mass, in order to search for variations larger 
than those expected from statistical fluctuations. These 
are evaluated using a method similar to the PDG scale 
factor |lll. Il5|. The corrections to the momentum scale 



TABLE I. Summary of the results from the fits to data for 
the D" — >■ A'^7r+ and D'^ — ^ K~n^Tr~n^ channels (statistical 
uncertainties only). S/B is the ratio of the convolved signal 
PDF to the background PDF at the given Am and u is the 
number of degrees of freedom. 



Parameter 


D° -^ Ktv 


D° -> Kmvn 


Number of signal events 


138 539 ±109 


174 286 ±150 


r[keV] 


83.5 ±1.7 


83.2 ±1.5 


scale factor, (1 + e) 


1.06 ±0.01 


1.08 ±0.01 


Amo[keV] 


145 425.6 ±0.6 


145 426.6 ±0.5 


S/B at peak 
(Am = 0.14542 [GeV]) 


2700 


1130 


S/B at tail 
(Am = 0.1554 [GeV]) 


0.8 


0.3 


2 / 


574/535 


556/535 



and dE/dx loss in detector material are varied to ac- 
count for the uncertainty on the Kg mass. We estimate 
uncertainties related to our choice of the radius assumed 
in the Blatt-Weisskopf form factor, and vary the reso- 
lution shapes according to their uncertainties and cor- 
relations. We vary the end point used in the fit, virhich 
affects whether events are assigned to the signal or back- 
ground component. This variation allows us to evaluate 
a systematic uncertainty associated with the background 
parametrization; within this systematic uncertainty, the 
residual plots shown in Fig. [1] are consistent with being 
entirely fiat. Additionally, we vary the description of the 
background distribution near threshold. Finally, we as- 
sess systematic uncertainties associated with the interval 
width and possible biases observed in fits to MC vali- 
dation samples with and without radiative effects. All 
these uncertainties are estimated independently for the 
D° — >■ K~TT^ and D'^ — >■ K~Tr'^TT~TT~^ modes, as dis- 



cussed in detail in Ref. [ll| and summarized in Table llll 

The largest systematic uncertainty arises from an ob- 
served sinusoidal dependence for Atoq on (f>. Varia- 
tions with the same signs and phases are seen for the 



reconstructed D'^ mass in both D'^ 



i?o 



tigation reveals that this dependence is most probably 
due to small sinusoidal variations in the magnetic field 
map used in track reconstruction [11| . The important 
aspect for this analysis is that the average value is un- 
biased by the variation in 0, which we verified using the 
reconstructed Kg mass value. The width does not dis- 
play a (p dependence, but each mode is assigned a small 
uncertainty because some deviations from uniformity are 
observed. The lack of a systematic variation of T with 
respect to (j) is notable because Atoq shows a clear depen- 
dence such that the results from the i)° — > A'^7r+ and 
D^ — >■ K~Tr~^Tr~TT~^ samples are highly correlated and 
shift together. We fit the Ar?io values with a sinusoidal 



function and take half of the amplitude as the estimate 
of the uncertainty. 

The results for the two independent 13° de- 
cay modes agree within their uncertainties. The 
dominant systematic uncertainty on the RBW pole 
position comes from the variation in (f) (1.5-1.9 
keV). For the decay mode D" — > K^n^ 

we find r = (83.4 ± 1.7 ± 1.5) keV and Ajtiq = 
(145 425.6 ± 0.6 ± 1.8) keV, while for the decay mode 
D^ -^ K-TT+'K-TT+ we find F = (83.2 ± 1.5 ± 2.6) keV 
and Amo = (145 426.6 ± 0.5 ± 2.0) keV. Account- 
ing for correlations, we obtain the combined mea- 
surement values F = (83.3 ± 1.3 ± 1.4) keV and 
Amo = (145 425.8 ± 0.5 ± 1.8) kcV. 

Using the relationship between the width and the 



coupling constant 11|, we can determine the experi- 
mental value of gD*+D°Tr+ ■ Using F and the masses 
from Ref. [l5| we determine the experimental coupling 
9'd*+d'>-k+ ~ 16.92 ± 0.13 ± 0.14, where we have ignored 
the electromagnetic contribution from D*+ — ?► D~^j. 
The universal coupling is directly related to go'Dv by 
9 = gD*+D°-K+f-Kl (2-y/m£)om£)*+). This parametriza- 
tion is different from that used by CLEO [8|; it is cho- 
sen to match a common convention in the context of 
chiral perturbation theory, as in Refs. [J, |l7|. With 
this relation and /^ = 130.41 MeV, we find ^'^'^p = 
0.570 ±0.004 ±0.005. 

Di Pierro and Eichten [18J present results in terms of 
i?, the ratio of the width of a given state to the univer- 
sal coupling constant. At the time of their publication, 
g = 0.82 ± 0.09 was consistent with the values from all of 
the modes in Ref. [il|. In 2010, BaBAR published much 
more precise mass and width results for the -Di(2460)° 

Using these values, our mea- 
we calculate 



and 1)2(2460)° mesons [li 



18| 



surement of F, and the ratios from Ref. 
new values for the coupling constant g. Table Hill shows 
the updated results. We estimate the uncertainty on g 
assuming err ^ T- The updated widths reveal signif- 
icant differences among the extracted values of g. The 
order of magnitude increase in precision of the D*"*" width 
measurement compared to previous studies confirms the 
observed inconsistency between the measured D*"*" width 
and the chiral quark model calculation by Di Pierro and 
Eichten ^. 
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TABLE II. Summary of systematic uncertainties with correlation p between the D" — > K vr^ and Z?" — >■ K 
The K~-K^ and K~-k^-k~-k^ invariant masses are denoted by m (D^cco)- 



Source 



a,j,.(r) [keV] 
Kit Krnvn 



Ktt 



:(Amo)[keV] 
Krnvn 



Disjoint D*+ momentum bin variation 

Disjoint m (Dreco) bin variation 

Disjoint azimuthal variation 

Magnetic field and material model 

Blatt-Weisskopf radius 

Variation of resolution shape parameters 

Am fit range 

Background shape near threshold 

Interval width for fit 

Bias from validation 

Radiative effects 



0.88 


0.98 


0.47 


0.24 


0.20 


0.28 


0.00 


1.53 


0.56 


0.04 


0.00 


0.22 


0.62 


0.92 


-0.04 


1.65 


1.81 


0.84 


0.29 


0.18 


0.98 


0.75 


0.81 


0.99 


0.04 


0.04 


0.99 


0.00 


0.00 


1.00 


0.41 


0.37 


0.00 


0.17 


0.16 


0.00 


0.83 


0.38 


-0.42 


0.08 


0.04 


0.35 


O.IO 


0.33 


1.00 


0.00 


0.00 


0.00 


0.00 


0.05 


0.99 


0.00 


0.00 


0.00 


0.00 


1.50 


0.00 


0.00 


0.00 


0.00 


0.25 


0.11 


0.00 


0.00 


0.00 


0.00 



Total 



1.5 



1.8 



2.0 



TT modes. 



TABLE III. Updated coupling constant values using the latest 
width measurements. Ratios are taken from Ref. Ilal. 



State 



Width (r) 



R = T/g^ 
(model) 



9 



D*(2010)+ 83.3±1.3±1.4keV 
Di(2420)° 31.4 ± 0.5 ± 1.3 MeV 
D2(2460)° 50.5 ±0.6 ±0.7 MeV 



143 keV 0.76 ±0.01 
16 MeV 1.40 ±0.03 
38 MeV 1.15 ±0.01 
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